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For the original hybrid inflation as well as the supersymmetric F-term and D-term hybrid mod- 
els, we calculate the level of non-gaussianities and the power spectrum of curvature perturbations 
generated during the waterfall, taking into account the contribution of entropic modes. We focus on 
the regime of mild waterfall, in which inflation continues for more than about 60 e-folds A'^ during 
the waterfall. We find that the associated /nl parameter goes typically from /nl — — l/Acxit in 
the regime with N ^ 60, where Aoxit is the number of e-folds between the time of Hubble exit 
of a pivot scale and the end of inflation, down to /nl ~ —0.3 when N > 60, i.e. much smaller 
in magnitude than the current bound from Planck. Considering only the adiabatic perturbations, 
the power spectrum is red, with a spectral index n^ — 1 ~ A/Na^it, in the case N ';$> 60, whereas 
in the case A > 60, it increases up to unity. Including the contribution of entropic modes does 
not change the observable predictions in the flrst case. However, in the second case, they are a 
relevant source for the power spectrum of curvature perturbations, of which the amplitude increases 
by several orders of magnitudes and can lead to black hole formation. We conclude that due to the 
important contribution of entropic modes, the parameter space leading to a mild waterfall phase is 
excluded by CMB observations for all the considered models. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

In the standard cosmological scenario, the large scale 
structures of the Universe are seeded by the quantum 
fluctuations of one or more scalar fields during a primor- 
dial phase of exponentially accelerated expansion. Due 
to this phase of inflation, the quantum field fluctuations 
becomes classical and are stretched outside the Hubble 
radius. They lead to a nearly scale invariant power spec- 
trum of curvature perturbations V(, whose amplitude 
and spectral index are constrained by CMB observations. 
The recent results from the Planck experiment [l| give 
the best bounds with In 10^° A^ = 3.089 ± 0.027 and 
n^ = 0.9603 ± 0.0073 with 68% C.L., where A^ is the 
value of Vi^ik^) at the scale fc* — O.OSMpc""^. Models 
of inflation predicting rig > 1 are now ruled out at more 
than the 5cr level [2]. 

This is commonly thought to be the case of original 
hybrid inflation Q. In this two-field model, infiation is 
realized in the false vacuum along a nearly flat valley of 
the held potential. It ends with a nearly instantaneous 
waterfall phase, triggered when the auxiliary field devel- 
ops a tachyonic instability. Hybrid models are well moti- 
vated from the point of view of particle physics because 
they can be embedded in frameworks like supersymme- 
try H-H, supergravity [lOl [ll| and Grand Unified Theo- 
ries |12h14J . Other realizations of hybrid inflation in the 
context of string theories |15l - [20| , and extra-dimensional 
theories [21| have also been proposed. 
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In addition to a blue tilted power spectrum, the orig- 
inal hybrid model was thought to suffer from a problem 
of fine-tuning of initial field values [23, l23|. Moreover, 
because a Z2 symmetry is broken at the end of inflation, 
it leads to the formation of domain walls with dramatic 
consequences for cosmology. 

However, recent developments have shown that there 
exists a regime in the parameter space of the original 
hybrid model for which the flnal waterfall phase is suf- 
flciently mild for inflation to continue for more than 60 
e-folds N after the flelds have crossed the critical insta- 
bility point [2J, I23]. In this scenario, topological defects 
are stretched outside the observable Universe, thus es- 
caping observation, and the power spectrum of adiabatic 
perturbations is red, possibly in agreement with CMB 
observations. Moreover, the initial fleld values leading to 
inflation have been shown to occupy an important pro- 
portion of the field space in a large part of the poten- 
tial parameter space |26l - [29| , thus solving the fine-tuning 
problem of initial conditions. Similar conclusions have 
been obtained for the most well-known supersymmetric 
realizations of the hybrid model, the so-called F-term and 
D-term models ^] . In the standard regime of a fast wa- 
terfall, the contribution of the cosmic strings formed at 
the critical instability point must be considered for these 
models. The F-term model was nevertheless found to be 
in tension with WMAP data whereas the D-term model 
was strongly disfavoured [3l|, [S^] ■ With the Planck re- 
sults, the degeneracy between the spectral index and the 
string tension, that tends to favor larger values for the 
spectral index, is strongly reduced and both models ap- 
pear to be ruled out in the standard regime [33|. 

The revival of the original hybrid scenario would 
be nevertheless of short duration if the mild water- 



fall scenario leads generically to a large level of non- 
gaussianities, with a local /nl parameter outside the 
recent Planck bound /nl = 2.7 ± 5.8 (68% C.L.) [H, 
or if the power spectrum of curvature perturbations is 
strongly affected by the entropic modes along the water- 
fall trajectories. 

The aim of this paper is to evaluate the power spec- 
trum of curvature perturbations and the contribution of 
entropic modes, as well as the level of non-gaussianities 
produced in the mild waterfall regime, for the original 
hybrid scenario as well as for the supersymmetric F-term 
and D-term models. Our calculation therefore differs 
from Refs. p35l - l44| since they consider only the scenario 
of a waterfall phase lasting no more than a few e-folds. 
A similar calculation was performed in Ref. ^a] for the 
original model, and some values for the parameters were 
found to lead to a large level of local /nl , but the focus 
there was on super-planckian field evolution in a regime 
where the potential is dominated by the separable terms. 
In this paper, we use an unified parametrization of the 
potential for the three considered models (F-term, D- 
term and original hybrid inflation) and restrict the anal- 
ysis to sub-planckian field values, such that supergravity 
corrections to the potential can be neglected. In this 
regime, the potential is not separable. 

We use the 5N formalism to calculate the local /nl 
parameter, both analytically and numerically, as well as 
the amplitude and spectral index of the power spectrum 
of curvature perturbations. In addition, because the SN 
formalism does not allow to study separately the time 
evolution of adiabatic and entropic perturbation modes, 
and since it cannot capture possible features in the power 
spectrum, we have also integrated numerically the lin- 
ear multi-field perturbations and derive the exact power 
spectrum of curvature perturbations. 

Denoting by iVexit the number of e-folds between the 
time of Hubble exit of a pivot scale and the end of in- 
flation, we flnd that the associated /nl parameter goes 
typically from /nl — —l/A^cxit, when the waterfall lasts 

for iV > 60 e-folds, down to /nl 0.3 when N > 60. 

In all cases the magnitude of /nl does not exceed the 
bounds of Planck. Considering only the adiabatic per- 
turbations, the power spectrum is red, with a spectral 
index Ug — 1 — 4/iVoxit, in the case A^ ^ 60, whereas in 
the case A > 60, it increases up to unity, which appar- 
ently suggests the presence of a parametric region, where 
the spectral index is in accordance with the observation 
by Planck [2^ |30|- When including the contribution 
of entropic modes, we flnd that the predictions do not 
change in the flrst case (A ^ 60). However, in the sec- 
ond case (A > 60), entropic modes are a sizeable source 
for the power spectrum of curvature perturbations and 
the spectral index first takes lower values before eventu- 
ally increasing up to unity. However, the amplitude of the 
power spectrum of curvature perturbations is enhanced 
by several orders of magnitudes and can even reach the 
level of black hole formation. We conclude that due to 
the important contribution of entropic modes, the param- 



eter space leading to a mild waterfall phase is excluded 
by CMB observations for all the considered models. 

This paper is organised as follows: in Sec.|ll]we give the 
exact background multi-field dynamics. Then are intro- 
duced the 5N formalism (Sec. IIII[) and the linear theory 
of multi-field perturbations fSec. lIVp . In Sec.|V]the con- 
sidered hybrid models are briefly described and the uni- 
fled parametrization of the potential is given. Sec. IVII is 
dedicated to the slow-roll dynamics during the mild wa- 
terfall phase. In Sees. IVIII and IVIlH we evaluate respec- 
tively the level of non-gaussianities and the power spec- 
trum of curvature perturbations, and compare the ana- 
lytical and numerical results. In the conclusion fSec.lIXp. 
we discuss the impact of our results on the constraints 
of hybrid models and envisage interesting perspectives to 
this work. 



II. MULTI-FIELD BACKGROUND DYNAMICS 

Assuming that the Universe was filled with n nearly 
homogeneous real scalar fields 4>i=i^2....n, the background 
dynamics is given by the Friedmann-Lemaitre equations 
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where H is the Hubble expansion rate, Mp\ = mpi/v87r 
is the reduced Planck mass, V^(0i=i,...,n) is the field po- 
tential, and where a dot denotes the derivative with re- 
spect to the cosmic time t. Then one can define a, the 
so-called adiabatic field [461] , that describes the collective 
evolution of all the fields along the classical trajectory, 
and the velocity field 



nS' 



(4) 



The equation of motion of the adiabatic field is given by 
'6 + 3Hcr + T4 = , (5) 

where 

dV 



2 — 1 



(6) 



with Uj being the components of a unit vector along the 
field trajectory Ui = (^ij a. 



III. THE 5N FORMALISM 

The 5N formalism, based on the separate universe 
approximation (and other assumptions discussed e.g. in 
Ref. [43]), states that the curvature perturbation C{x,t) 
on a spatial hypersurface of uniform energy density is 
given by the difference between the number of e-folds 
realized from an initially flat hypersurface N(t, x) = 
ln[a{t)/a{ti)], where a{x,t) is the local scale factor, and 
the unperturbed number of e-folds No{t) = \n[a{t)/a{ti)]. 
If we label respectively by i and f the initial and final hy- 
persurface, one has 



C = SN^ = Nit,x)-Noit) 



(7) 



Our initial hypersurface is chosen at the time i* corre- 
sponding to the Hubble exit of the observable pivot scale 
A:* — 0.05 Mpc~ (in the rest of the paper, star subscripts 
indicate quantities evaluated at t^). The final hypersur- 
face must be of uniform energy density. In this paper we 
are interested in the curvature perturbations at the end 
of the slow-roll regime, when one of the slow-roll param- 
eters reaches unity. 

The field perturbations are close to Gaussian and have 
a very small amplitude, so that the observed curvature 
perturbations are given to good accuracy by 



n 1 ^^ 






where we have used the notation 
dSN^ 



N, = ^,Nii 
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The amplitude of the reduced bispectrum is given by [48| 

(10) 
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One can also calculate the power spectrum amplitude 
and spectral tilt as 
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where ei = —H/H^ is the first slow-roll parameter. Prac- 
tically, instead of a final surface of uniform density, we 
have considered a final surface of constant slow-roll pa- 
rameter value, the e-fold difference between these two 
surfaces being negligible compared to the e-folds differ- 
ences generated by varying the initial field values. The 
validity of this assumption has been checked numerically. 



We furthermore need to introduce the number of e- 
folds A^* in the sense of a reparameterization of time, 
dN* = Hdt, implying that during de Sitter inflation, 
where a{t) — exp{Ht), N* — Ht. We use the subscript k 
to indicate its value when a given scale k exits the hori- 
zon, * when this scale is the pivot scale and the subscript 
'end' to indicate its value at a uniform density surface 
at the end of inflation. As stated above, the latter can 
be practically obtained by checking for the violation of 
the slow-roll conditions, c/. also the discussion in Sec- 
tion IVTlBj Finally we deflne the function 



N ^ Nl^^ 



iV* 



(13) 



giving the number of e-folds up to the end of inflation. 



IV. THE LINEAR THEORY OF MULTI-FIELD 
PERTURBATIONS 

A. Perturbed equations 

The perturbed metric can be written as 

ds^ ^ a^ [-(1 + 2$)d772 + (1 - 2*)%da;Ma;^] , (14) 

where $ and Vt are the Bardeen potentials and rj is the 
conformal time, which is related to the cosmic time t as 
dt — adrj. In the longitudinal gauge, the spatial non- 
diagonal Einstein equations perturbed at first order lead 
to $ = ^, such that the (0,0), (0,i) and {i,i) equations 
read 

-3-H($' + 'H$)+V2$ 
-^E{'^'M-€'^ + a'^^S<^?j , (15) 
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where H = a' /a and Scjji is the perturbation of the scalar 
field (j)i, and where a prime denotes the derivative with 
respect to the conformal time. On the other hand, the n 
perturbed Klein-Gordon equations read 



^" + 2HS(I)', - V'^ScI), + Y, a^Hj 
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The field perturbations are coupled to each other through 
the cross derivatives of the potential and the Bardeen 
potential. By adding Eq. (jisp to Eq. (fT7|). and by us- 
ing Eq. (J16p . one obtains the evolution equation for the 



Bardeen potential, 



4-H^)$- V^* 
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We want to obtain the curvature perturbation C, defined 
as 

H 



C=<^- 



H' -H^ 



($' + -H$) 



(20) 



The background dynamics imphes that T-L' — H^ = 
— 47r(T'^/mpj. By using Eq. (1161) . the comoving curvature 
can thus be rewritten 



c = $ 



T-l 
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4=1 



From the background and the perturbed Einstein equa- 
tions, one can show that C, evolves according to [431 
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UiUj has been 



where the orthogonal projector _Ly = Id - 
introduced. For a single field model, the second term 
vanishes and one recovers the one-field evolution of C,. 
For the multi-field case, one sees that entropy perturba- 
tions orthogonal to the field trajectory can be a source 
of curvature perturbations, even after Hubble exit. 



B. Numerical integration 

For the numerical integration of the perturbations, we 
refer to Ref. ^49|| and give here only the guidelines for 
the calculation of the exact power spectrum of curvature 
perturbations in a multi-field scenario. It is convenient 
to use the number of e-folds as the time variable. Some 
equations are redundant, and one can for instance use the 
Bardeen potential expressed directly in terms of the field 
perturbations 5(j)i and their derivatives. After expanding 
in Fourier modes, Eq. ([18]) reads [4^ 
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Here, S(f)i 
isotropy. 



hik, rj) and we may use k = |k| because of 



The quantization of the field perturbations in the limit 
k ^ aH provides initial conditions for the 5(j)i. For the 
field operator, we take 



(50i(?7,x) 



(fik 
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a,(k)e-*''-''(5(/),(fc,r/)+h.c. , (25) 



where h.c. stands for hermitian conjugation and 
[ai(k),a](k')] — (27r)^5jj(5(k— k'). The normalized quan- 
tum modes are defined by 



Vi.k{ri) = aS4>i{k,r]) . 



(26) 



Neglecting the mass terms, they obey the equation w"^ -|- 
k'^Vi,k — 0, and in the regime k ^ aH , one has 

lim Vk4r]) = fce-^'^^"""') . (27) 

In terms of the field perturbations, the initial conditions 
(denoted by the subscript i.e. to avoid confusion with 
previous notation) therefore read, up to a phase factor. 
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It is not convenient to integrate the perturbations from 
the onset of inflation, since the total number of e-folds 
can be much larger than A*. In order to avoid the time 
consuming numerical integration of sub-Hubble modes 
behaving like plane waves, it is convenient to start to 
integrate the perturbations later, when the condition 



k 



H(ni.e.) 



= Cfc » 1 



(30) 



is satisfied, where Ck is a constant characterizing the de- 
coupling limit. To summarize, the numerical integration 
of multi-field perturbations can be divided in four steps: 

1. The background dynamics is integrated until the 
end of inflation, such that Nl^^ and Nl^^ — Nl are 
obtained. 

2. The background dynamics is integrated again, until 
A'j*^, is reached. Initial conditions for the perturba- 
tions are flxed at this time. 

3. For each comoving mode fc, the background and 
the perturbation dynamics are integrated simulta- 
neously from A^*^ to N^^^^. 

4. Determination of the scalar power spectrum, 
•Pj-(fc) = l/{2iT'^)Y.i\Ci?^ where i = l...n stands 
for the n independent initial conditions for the held 
perturbations 50i, Q being the induced contribu- 
tions to the curvature perturbation ^. 



HYBRID MODELS 



A. Original version 



The original hybrid model of inflation was first pro- 
posed by Linde ^3,]. Its potential reads 



T/(0,V) = A 
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The field is the inflaton, ?/; is an auxiliary transverse 
field and M, /x, (pc are three parameters of mass dimen- 
sion. Inflation is assumed to be realized in the false vac- 
uum [5] along the valley {ip) = 0. In the usual descrip- 
tion, inflation ends when the transverse fleld develops a 
Higgs-type tachyonic instability soon after the inflaton 
reaches the critical value 0c. From this point, the classi- 
cal system is assumed to evolve quickly toward one of its 
true minima (0) = 0, ("0) = ±M, whereas in a realistic 
scenario one ex pec ts the instability to trigger a tachyonic 
preheating era |50l - [55j . 



B. Supersymmetric F-term model 

Supersymmetric i^-term inflation has been proposed 
in Refs. [J, |5| and has subsequently been discussed ex- 
tensively in the literature. The underlaying details for 
our present analysis are given in Ref. [30|, and we only 
present the main assumptions and results for the poten- 
tials here. The superpotential is given by 



W = nS{HH - m2 



(32) 



where the superfleld S* is a gauge singlet and the super- 
flelds H (H) transform in the (anti-)fundamental rep- 
resentation of SU(A/'). This implies a tree- level scalar 
potential 
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where now 5, H and H are complex scalar flelds. The 
i'-term gives a large mass to the field combination 
{l/\/2){H — H), such that we do not need to consider 
this field direction for the dynamics of fluctuations dur- 
ing inflation. The waterfall field can be identified as 
ip = {l/y/2){H + H)) and the inflaton fleld as = y/2\S\. 
In terms of these degrees of freedom, the tree-level po- 
tential is 
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In the potential valley where > 0c = \/2m and ("0) = 0, 
the potential energy A = K^m^ spontaneously breaks su- 



persymmetry. At one-loop order, this leads to the follow- 
ing correction to the potential: 
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is a renormalization scale. 
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(36) 



where Q is a renormalization scale. Notice that the 
derivatives of this potential with respect to that are 
phenomenologically relevant for inflation are independent 

ofg. 

We are concerned in this work with the dynamics near 
the critical point where w 0c and ^ sa 0. As discussed 
in Ref. 30] , treating the one- loop potential to linear order 
is a good approximation in that regime. Of relevance is 
the first derivative of the one-loop potential. 
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In order to realize a substantial amount of e-folds below 
the critical point, the relation k <^ 02 /Mp must hold (30| . 
It can then easily be derived that the second derivatives 
of the one-loop potential may be neglected and that they 
do not yield a phenomenologically relevant contribution 
to the spectral tilt of the observable power spectrum 130| . 



C. Supersymmectric D-term model 



Here, the superpotential is 



W = kSHH , 



(38) 



where S is again a singlet and H and H transform ac- 
cording to the one-dimensional representation of a U(l) 
gauge group. Spontaneous symmetry breaking is induced 
by the D-term 
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where mpi is the Fayet-Iliopoulos term. The infiaton 
field is again = a/2|S'|, whereas now the waterfall field 
is given by -0 = V2|-ff|- The critical point, where ip 
becomes tachyonically unstable is here 
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In terms of these various fields and variables, the tree- 
level potential can be expressed as 
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The one-loop potential is readily obtained from its ex- 
pression in the F-term case ([55)) when setting Af — 1. 



VI. INFLATION ALONG WATERFALL 
TRAJECTORIES 



D. Unified parametrization 

In this paper, we study the original hybrid model as 
well as the supersymmetric F-term and D-term vari- 
ants in a unified approach, by considering the following 
parametrization of the two-field potential. 



V{^,^)=A 
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where M, (j)c are respectively the position of the global 
minima and of the critical point of instability along the 
valley. In the case of the original model, one has p = 2 
whereas the dynamics for the F-term and D-term models 
near the instability point is well described when setting 
P = 1 [sol ■ The relations between the parameters of this 
potential and the model parameters for F-term and D- 
term inflation are given in TABLE HI 
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TABLE L Parameters to be substituted into the potential (I42p 
in order to obtain the F- and D-term models close to the 
critical point. 

For the purpose of deriving the phenomenological con- 
sequences of hybrid inflation, it is useful to note the 
derivatives 
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In particular, the first derivatives enter the slow-roll 
equations of motion, 

dV ■ dV o V 

d(j> dip 3Mpj 

We also make use of the standard definition for the slow- 
roh parameters tjxy = M^^[d^V/{dXdY)]/V, where X 
and Y can be either of the canonically normalized fields 

(/> and ip. 



We are interested in the field dynamics during the wa- 
terfall regime, i.e. in the times after the field trajec- 
tories cross the critical instability point 0c- Tachyonic 
preheating is triggered when the exponentially growing 
long- wavelength perturbations of the auxiliary field ip be- 
come non-linear. This only occurs when the tachyonic 
auxiliary field mass is larger than the Hubble expan- 
sion rate, m'^, > H^ ~ V/{3M^{), or equivalently when 
—V^^> <^ 1- In the commonly studied large- coupling sce- 
narios of hybrid infiation, this condition is satisfied after 
no more than a few e-folds after entering the waterfall 
phase. However, there exists a large region in the pa- 
rameter space (the small coupling limit for the F-term 
and D-term models) for which the waterfall phase lasts 
for more than 60 e-folds. We refer to this situation as the 
mild waterfall. The field evolution in this regime has been 
studied for the original hybrid model in Refs. [24.I25J. and 
for the F-term and D-term models in Ref. [30|. In this 
Section, we re-derive the mild waterfall dynamics for the 
unified potential of Eq. (H^ . For the purpose of an ana- 
lytic description, we follow Ref. \2^ and divide the field 
evolution below the critical point into three phases: 



• Phase 0: when the second term of Eq. (I43bp and 
the first term of Eq. (I43ap dominate. 



Phase 1: when the first term of Eq. (|43bp and the 
first term of Eq. (|43ap dominate. 



• Phase 2: when the first term of Eq. (j43bp and the 
second term of Eq. (|43ap dominate. 

Classically, the phase does not last more than about one 
e-fold [25;] . Moreover, in a realistic scenario, the quantum 
diffusion of the auxiliary field dominates over the classical 
dynamics very close to the critical instability point. For 
these reasons we do not consider the phase but only 
the phases 1 and 2. 

Following the derivation of Ref. [25|, it is convenient 
to parameterize the fields as 
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where ipQ is the initial condition for the auxiliary field 
at the critical point of instability. During the slow-roll 
waterfall regime, one has f < and |^| ^ 1, which is 
consistently verified by the explicit solutions. 



A. Phase 1 

Solving the slow-roll equations, one finds that the field 
trajectories during the phase 1 follow the relation 
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as long as the temporal minimum (corresponding to the 
ellipse where dV/dip vanishes) is not reached. In this 
case, phase 1 connects to phase 2 at the point (^2,X2), 
with 



X2 = ln 






(47) 



One obtains the number of e-folds realized in phase 1 by 
integrating 



diV* 






(48) 



Imposing that iV* = at the critical point, one finds 



e = - 



pK€-' .rt 



fiP 



-N' 



(49) 



The temporal minimum for the field ^p is located on the 
trajectories 



C 



2p2x 



%e' 



2M2 



(50) 



(51) 



It is reached during phase 1 if the condition 

^^ 16/iPM2 
is satisfied. One therefore gets 



_ Vp^rM ^^^ 

f, = ^ 2/iP/2 ^^2 I"'^ X2 > ^Q^p^p .^^ 

^ P^c r ^ prc+' 

^ tor Y2 < ^ 

8fiP 16AiPAf2 

The number of e-folds iVi realized in phase 1 then follows 
as 



^^~~i ^°'' ^2 > -.^ „,,o 

iVr = 2v^^|-M2, 16m^A 2 _ 

P0P+2 ^ ^ 



8M2^ 



for X2 < 



IGfiPM 



,T2 



B. Phase 2 



The slow-roll equations in phase 2 yield the trajectory 
^ pAf20P-2 



C^ = ^2^ 



8/iP 



.2(x-X2) 



1 



(54) 



Once again, we have to distinguish between the possibil- 
ity where this trajectory reaches the temporal minimum 
at some point before the end of inflation, and the possibil- 
ity where the temporal minumum is not reached during 



inflation. In phase 2, the temporal minimum is reached 
when 



e = e 



2T.M. 



M1^JM1 + .2_P^ .55^ 

202 V40^ ^^ S,.P- (^^) 



On the other hand, inflation ends before reaching the 
temporal minimum when 



VipTp 



8M^ 
M2 



-1 



(56) 



In the opposite case, where the temporal minimum is 
reached before the end of inflation, the slow-roll condi- 
tions are violated when 



4Mp2j^2 



1. 



■'W- ^2M2 

By using Eq. ((50)) . one thus finds that 



(57) 



->cnd 



8M^, 
8M| 



for l^cndl > |6t.m. 
for |<^ond| < |6t.m. 



(58) 



During phase 2 and before reaching the temporal mini- 
mum, assuming that ^ ^ 1 and X2 > 1/2, the slow-roll 
equations in c-fold time can be solved exactly and one 
finds [2i 



m') 



I pAP {c' ~c)f{N*)~c' -, 
2<j,r^fiP (c' - c)f{m) +c' + , 



(59) 



where c = ^^2/2, c' = y/p — 1/4 and 



/(iV*) = exp 



8V2c 



^Af2(iV*, 



iVi) 



^Jpp,PAP 



(60) 



A good approximation can be obtained by considering 
the limit |^| ^ |^2|, where one can obtain 



1 1 



8M| 
A/2 



(iV* 



^end) 



(61) 



The number of e-folds realised along the temporal min- 
imum is obtained by integrating 



diV* 



This gives 



^c„d 



/v* 



8M2,^ 



In 



8Af2j 



•^cnd 



2T.M. 



(62) 



(63) 



where iVj^p y^ is the value of the parameter TV* when 
^ = C2T.M.- If we understand the model specified by 
the potential (|42)) as an effective theory valid below the 
Planck scale, we should impose that 0^ <C Afpi and con- 
sequently, the number of e-folds inflation continues for 
after reaching the temporal minimum is very low. 
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VII. NON GAUSSIANITIES FROM THE MILD 
WATERFALL PHASE 

A. Analytical results 

In this Section, we calculate the level of non- 
gaussianities using the SN formalism. We denote the 
time where the transition between the phases 1 and 2 
occurs by ti^2 and the time of Hubble exit of the pivot 
scale fc* by t*. We need to distinguish the case when 
the pivot scale leaves the Hubble radius in phase 2 (case 
1: long waterfall phase with A^ond ^ ^0, and t* > ^1,2) 
from the case where horizon exit occurs in phase 1 (case 
2: moderately long waterfall phase, with Nj:^^ > 60 and 
t* < ^1.2)- 



1. Case 1: Hubble exit in phase 2 

The first step for calculating the local /nl parameter 
with the 5N formalism is to derive the number of e-folds 
until the end of inflation starting from an arbitrary ini- 
tial point (^j, Xi) in field space. In order to compare with 
the observational bounds, we then use this result to eval- 
uate the quantities N^^, N,^, N^(f,(f,, N^^^ and N^^^ for 
the pivot scale /c* , at which the unperturbed background 
fields take the values (^»,x*)- 

We first consider the case where inflation ends before 
the temporal minimum is reached. Integrating the slow- 
roll equation of motion in phase 2 give the trajectories 



^^^^^^e'x^24>i{e~m+ro^ 



2p2x, 



(64) 



(This equation also determines Xcnd when replacing ^ -^ 
^ond and X ~^ Xond-) Above relation can be used in com- 
bination with the slow-roll equation for </> to derive the 
number of e-folds that elapse while the fields evolve from 

(?i,Xi) to (Ccnd,Xcnd), 



^end 



iv; 



AP_ 



d^ 



e 






(65) 



Defining C = -Q + tpoc^^^ / {2(j)l) , one finds 

(1^2 



NL<x -Nl = - 


IVl 




SM^.VC 




X 


arctan — ^= 

\VcJ 


- arctan ^^ 


if C > 0, and 




N* .- N* = 

-"end ^'j 


A'P 




8M^,VC 




X 


, ,/ecnd\ 

arctanh --= 

I \VcJ 


— arctanh —^ 



(66) 



(67) 



if C < 0. 

For applying the (5iV formalism, we need to evaluate 
the e-fold derivatives with respect to 0i and V'i , evaluated 
in the limit where ^i — > ^ and Xi ^ X^ where ^ and x 
belong to the unperturbed waterfall trajectory, and for 
which the relations 



2e2x 



V-oe 



and 



2<^cC' 



M2 



8M2j 



end 



m 



M2 



SA-I^i 



(68) 



(69) 



are satisfied. Note that in this limit, one gets C — >■ and 
one can expand Eqs. (|66p and (1571) in a Taylor series to 
obtain 



TV* 



N\ 



a£_ 
Mil, 



1 



1 



M 
3e? 






end 



(70) 



which is a consistent generalization of Eq. (|6T|) . We 
use this expression to calculate the e-fold derivatives 
N^if,, N^^, N^ij,ij,, N^^jf, and iV_^^ around the field configura- 
tion (^i, Xi)- For the purpose of calculating these deriva- 
tives, we can now relabel (Ci,Xi) ^ (C;X)- P-e- £, has 
now a different purpose than the integration variable in 
Eq. ([65)) ]. We thus obtain 



N^ = 



NP 



ce^Ki 



1 

e 



2 



2g 
3 

end 



3^, 



M' 



2f2 



240cMp^ie 



(71) 



and 



N^ = 



AP_ 
8M^ 



^ 



V- 



30?e 30?^, 



end 



M'^i) 



2« 



2402Mp^,e 



(72) 



The approximations above are valid when |^| < |^ond| 
and thus |C'^| ^ I Condi- ^s ^^ eventually replace f — >■ ^*, 
these relations are well satisfied for the present purposes. 
By using ^ = — -\/2(/>cC, we notice the useful relation 



N_, 



V2 



N_, 



(73) 



For the second derivatives, when keeping only the lead- 
ing terms, we find 



N 



A'P 



Ip-ip 



2mMi, 



N_, 



1 

e 



1 



e 



AP 



end 



ANP 



1 c-i 



24(/.2Mp^ie 



2f3 



240^Mp-,e 



(74) 



(75) 



and 



^,<t>il> ~ 



M^^k 



ble^Ml.e 



3\/2M2 



(76) 



One can now evaluate the /nl parameter at the pivot 
scale fc*. We first notice that 



N_, 



AN 



■TpTp 



3^2 



N 



(flip ■ 



(77) 



Using that the number of e-folds inflation lasts after the 
horizon exit of the pivot scale is given by iVcxit = ^cnd ~ 
Nl and evaluating above derivatives for (^, x) = {£,*,X*): 
we obtain 



/i 



5 /24M2ie, 

NL — 777 I - 



18 



M2 



HN 



exit 






<1 



(78) 
In the case of the F-term and D-terni models, this ex- 
pression reduces to 



/i 



NL 



3(iVexit + 1) 



(79) 



The level of non-gaussianities is therefore negative and 
very low, typically /nl ~ —0.03 for the F-term and D- 
term models. 

We finally consider the situation where the temporal 
minimum for the field ip is reached in phase 2 before the 
end of infiation. Then, we need to evaluate 



N = N 



T.M. 



N 



T.M. 



(80) 



The first term on the right hand side is the number of e- 
folds realised between t^ and the time where the temporal 
minimum is reached. The second term is the number of 
e-folds realised along the temporal minimum up to the 
end of inflation. 

Starting from an arbitrary initial field value in phase 
2, one finds the value of ^ at which the trajectory crosses 
the temporal minimum, 



S2T.M. 




i^t/im) 



(81) 



For evaluating NJ^-^- , one can take Eq. (f70| and replace 
Ccnd by ^2T.M.- For the derivatives, we again replace 
{£,i,Xi) -> (Ci x) and notice that 



d^ 



2T.M. 



< 



dC 



20? 



6 



(82) 



2T.M. 



This gives an additional term — (d^2T.M./d(/)i)/^x m i'^ 
the round brackets of Eq. ([7T|). that is negligible com- 
pared to the leading term. Next, using Eq. ([S5)) . one can 
calculate 



-"^T.M.. 



d6 



2T.M. 



8Mj^i6t.m. d^ 



M, 



2 0(0- 

pi 



(83) 



Since 1^1 < le 



iend 



Nj-}^-, provided 



cjJ^JiSM^i), we find that iV|ti.,^ « 
^ J , jjivjviucu. ^c <C Mpi, as one should require for the 
effective theory description to be valid. 

In a similar way, one finds that the leading terms in 
N^'^M.,^, ^t"a'/.,^^,.^t"^/.,0^ and iV^",^,.,^,^, as weh as 
in /nlj are not modified when the temporal minimum is 
reached, apart that one has to replace fond by ^t.m. • As 
a consequence, the level of non-gaussianities is reduced 
compared to the contributions from iVj*^' and cannot 
in any case increase up to an observable level. 



2. Case 2: Hubble exit in phase 1 

Now we consider the situation where the pivot scale 
exits the Hubble radius in the phase 1. We first consider 
the case where the temporal minimum is not reached. In 
order to obtain the spectra using the SN approach, we 
need to generalize the analysis of Section FVlAl where we 
have fixed the initial value of f to zero, to more general 
initial values fi, such that we can obtain the necessary 
e-fold derivatives. 

By integrating the slow-roll equations of motion, we 
obtain the trajectory 



e - e? 



4^P 



(84) 



Eq. (H51) gives the number of e-folds elapsing from the 
point {^i,Xi) until reaching (^217X2)1 where phase 1 ends. 



Nl 



Nl 



m 






plV 



AP-2 



(85) 



where ^2^ = --y/Cf +P</'? '^M^{x2 - Xi)/(4MP), Nf is the 
value of the parameter A^* at the initial point (fi, Xi)? -^2 
its value at (■^217X2), and X2 is again given by Eq. P7|) . 

The number of e-folds between the point (fiiXi) and 
the point where inflation ends in phase 2 due to the viola- 
tion of the slow-roll conditions is then given by Ni + N2, 
where A'2 are the e-folds between the onset of phase 2 
(^217X2) and the violation of slow roll. For the following 
purposes, we can again drop the index i on the initial 
field values, i.e. (C»,Xi) -^ (C,x) and ((/>i,V'i) ^ ('A,V')- 
The e-fold first derivatives then read 



N_, 



{N 



i,c- 



■N 



2,e; 7 



Vk 

We first evaluate Eq. (|57)) . The first term gives 
•P di2r M^ 



(86) 



(87) 



^^1.. 



/^^ 

r^Mi dx Mii.i 



'pi 
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and the second term yields 

p<j>P-'^M^ dN2 



N, 



2,X 



le^A^p d6 



(89) 



Using Eq. (|70|) with ^2i and X2i instead of ^i and Xi for 
the number of e- folds reahzed in phase 2, one gets 



diVa 



M2 



1 „ 2 , 2^; 



AP 



and thus one obtains 



^2,x = 



384e2^MpVP 12x2 



-A^i 



1,X: 



(90) 



(91) 



where the last expression is found by inserting the value 
of ^2 from Eq. (15^ . Since X2 > 0(1) [it is a large loga- 
rithm cf. Eq. (j47p ]. the dominant term in Eq. (|87p comes 
from the number of e- folds realised during phase 1. In 
the following, we will therefore neglect the term N2.X, 
such that 



N, 



M' 



8Mp^i6^ 



(92) 



We now turn to Eq. 

dJVi _ yf (e. 

dC pA/p>r' U2 

and for the second term, 

dA^2 d7V2 C M^S, 



For the first term, one finds 



dC de2. 6 24M2 ^1 



(93) 



(94) 



Using the appropriate expression from Eq. ([5^ for ^|, 
one sees that 



dA^i „ /6 A d7V2 



(95) 



Hence again, we may neglect the contribution coming 
from the number of e-folds in phase 2 and obtain from 
Eq. dMl) 



N_,, 



fi" 



pM^,<pr' 



(96) 



Moreover, since ip* *C V'2, we may use that A^^^ ^ N^jp. 

For the second derivatives, one can again show that 
the dominant terms come from the number of e-folds in 
phase 1. One therefore obtains 



N 



AP 



,tptp 



^Mi.i^r 



M2 



1 



1 

2x^ 



SMp^i^V' ^ 



N_, 



(97) 



TV 



</"/' 



M2$ 



8M2ei0c^ Cl'/'c 



^., 



and 



A^.. 



pMp>r'<^2 



C2 6 el 



pMpVr'^'/'? 



(98) 



(99) 



(100) 



Now we can evaluate the parameter /nl in the limit 
where the derivatives with respect to the field (j) as well 
as the contribution from the number of efolds in phase 2 
are negligible. One finds 

f2, 



/nl 



5 N^^jN^)^ 20Mp^i6 
'6 {N^f 



3Af2 



3M/iP/2 
For the F-term model, this implies 



(101) 



1/1 



5KVA^^/Tog2M, 



NLl 



6\/27i 



P'V^2<0.13VX^ (102) 



where the maximal negative value of /nl is obtained 
when the waterfall lasts just about 60 e-folds, i.e. when 
KvAf « m^/M^i [30|. Particle physics experiments im- 
pose a lower bound on tti, inducing the upper bound 
yX2 ^ 6. The negative amplitude of the /nl parame- 
ter is therefore never higher than about unity, which is 
below the Planck sensitivity. We have plotted in FIG. [T] 
the level of local non-gaussianities given by Eqs. (|10ip 
and ([78|) . as a function of k for the F-term model with 
different values of m. For the D-term model, the level of 
non-gaussianities is given by the same expression, with 
TOpi instead of tti. It is independent of the model param- 
eter g, apart logarithmically through the y'x2 factor. 

For the original model, requiring 60 e-folds along the 
waterfall imposes MiJ, > 40Mp[ |2J|. The level of non- 
gaussianities is independent of the parameters A and (j)c 
(apart through ^yx2), and its maximal negative value is 
therefore of about /nl — —0.3, as for the F-term and 
D-term models. In TABLE |ll] we give the level of non- 
gaussianities from the analytical approximations in the 
original model for various parameter sets, covering the 
qualitatively different regimes, and compare to the nu- 
merical results. Notice finally that the particular case 
where the temporal minimum is reached in phase 1 is 
not relevant, because that would imply that X2 ^ 1 and 
thus the quantum diffusion would still be dominating at 
the time of Hubble exit of the observable scales. 



B. Numerical analysis 

The analytical results of the previous Section are valid 
under some approximations, namely: i) the slow-roll ap- 
proximation, ii) a sharp transition between phase 1 and 
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phase 2, iii) the final hypersurface given by the condi- 
tion ^ = ^ond [according to Eq. ([55|] instead of a uniform 
density condition, iv) the Taylor expansion (|70|) of the 
number of e- folds in phase 2, v) some terms neglected in 
the e-folds derivatives. To check the validity and the ac- 
curacy of our results, we have implemented a numerical 
calculation that makes use of the ^A^ formalism. Practi- 
cally, we use the following algorithm: 

1. First, a numerical integration of the exact multi- 
field background dynamics, from the critical insta- 
bility point, where we set iV* = 0, until the end 
of infiation in order to obtain the total number of 
e-folds -/Vpnd of inflation along that part of the field 
trajectory. 

2. Second, an integration of the background dynamics 
from the critical point down to the time f * of Hub- 
ble exit of the pivot scale fc* = O.OSMpc^^. With 
the help of the first step, we know that the exit 
point is reached when the e-fold parameter equals 



*cxit ■ 



3. Numerical integration of the field dynamics from 
initial conditions on a 3 x 3 grid of values centered 
on (0* ,■)/'*) . Determination of the number of e-folds 
N to reach the final hypersurface. 

4. Numerical evaluation of the derivatives N^i and 
N^ij. The /nl parameter as well as the amplitude 
V^ and the spectral tilt Ug of the power spectrum 
of curvature perturbations can then be computed 
by using Eqs. ^, ^ and (fT2|). 

For the stability of the code, notice that the differences 
between the initial conditions need to be carefully chosen, 
sufficiently small for the numerical derivatives N^i and 
N^ij to be accurate, but at the same time sufficiently 
large for the differences between the values of N to be 
much larger than the integration steps (that cannot be 
lower than AA^ ~ 10~* without increasing unreasonably 
the integration computing time). 

Another numerical issue is related to the choice of the 
final hypersurface. It is particularly tricky to define it 
numerically because the variation of the false vacuum 
potential along the waterfall trajectories is so tiny that 
it cannot be resolved due to the limited numerical preci- 
sion. We have thus considered the following alternatives: 
i) by defining the final hypersurface as the end of the 
slow-roll regime, more precisely when one of the slow- 
roll parameters reaches unity, ii) by implementing the 
uniform energy condition after the end of inflation, when 
for instance pond = 0.99pinf (in this case the effects of the 
tachyonic preheating are considered to be negligible and 
the classical trajectories are assumed to be valid down to 
the final hypersurface) , iii) by defining a hypersurface of 
uniform potential energy density without including the 
constant term Vq = A. We did not find noticeable differ- 
ences between these possible methods. 



At the critical instability point, the classical dynam- 
ics of the auxiliary field is dominated by the quantum 
diffusion. Classical trajectories must therefore be seen 
as emerging from this quantum stochastic regime. They 
are valid once the classical vacuum expectation value of 

the auxiliary field fulfiUs ip ^^ \/{''Pqu)j where V'qu is the 
operator for the quantum field fiuctuations around the 
classical expectation value ip. For the numerical integra- 
tion of the field dynamics, we take for the initial condition 
of the auxiliary field at the critical instability point the 
value given by the process of quantum diffusion [2J, [30| 



V'O 



ipo 



-ipo 



2\/37r3/4(ln2)i/4 



for the F-term , 



(103) 



gmli 



8\/3^7r3/4(ln2)i/4 



for the D-term , (104) 



A/iM 



for the original model 



(105) 



967r3/2 
and only consider the trajectory as valid when ip '^ ipo- 

For the F-term model, we have plotted in Fig.[l]the lo- 
cal /nl parameter calculated numerically as a function k, 
for different values of the mass parameter m in the range 
10~''Mpi < M < O.lAfpi. For comparison, the analytical 
approximations derived in the previous Section are also 
presented. Numerical and analytical results agree well, 
with /nl values ranging typically from —0.03 down to 
—0.3, as long as m ^ O.OlAfpi. For larger values, in the 
regime i* < ti^2, our approximation given by Eq. (jlOip is 
not valid any more and the |/nl| parameter takes lower 
values. This is because the e-folds derivatives with re- 
spect to the field (p become an important contribution 
to the level of non-gaussianities. We can draw identi- 
cal conclusions for the D-term model, where the level of 
non-gaussianitiy is nearly independent of the coupling g. 

For the original model, the parameter space has two 
additional dimensions and it is more difficult to explore it 
entirely with our numerical method. However, we present 
in TABLE |ll] a comparison of numerical and analytical 
results for a few sets of parameters corresponding to the 
different regimes. This serves as a check that the numer- 
ical results are in agreement with the analytical approx- 
imations. 



VIII. POWER SPECTRUM OF CURVATURE 

PERTURBATIONS AND CONTRIBUTION OF 

ENTROPIC MODES 

A. Using the SN formalism 

The amplitude and spectral index of the power spec- 
trum of curvature perturbations at the end of inflation 
can be calculated by using the SN formalism, as ex- 
plained in Section [HII For the spectral index, in addition 
to the e-folds derivatives that have been calculated in the 
previous Section, one needs to know the fleld derivatives 
at the time of Hubble exit of the pivot scale. 
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Parameters 


Regime 


rnum 

Jnl 


rapp 
JNL 


nr^' 


nt"" 


Pr"(fc*) 


P^.^'-ik,) 


M = 0c = 0.01 A/pi, fi = lO'' Mpi, A = lO-''" M*i 


ii.2 < i* 


-0.030 


-0.029 


0.929 


0.930 


0.017 


0.019 


M = (?ic = 0.001 Mpi,/i = 10« Mpi, A = 10^^" M;Ji 


tl,2 < U 


-0.033 


-0.033 


0.921 


0.921 


1.1 X 10-^ 


1.3 X 10-'' 


M = 0.001 Mpi,<^c = 10-^/i = 10** Mpi, A = 10-3« M^i 


tl,2 < U 


-0.033 


-0.032 


0.921 


0.921 


0.0011 


0.012 


M = 0.01 Mpi, (^c = 0.1 Mpi, /i = 10^ Mpi, A = 10"^° M^^ 


tl,2 < tf < Itm 


-0.030 


-0.030 


0.929 


0.929 


1.8 X 10-* 


1.9 X 10-* 


M = 0.001 Mpi, 0c = 0.1 Mpi, ^l = 10'' Mpi, A = 10"^° M^^ 


il,2 < t» < tTM 


-0.033 


-0.032 


0.921 


0.921 


1.1 X 10-" 


1.3 X 10-" 


M = 0c = 0.01 A/pi, ^ = lO'' Mpi, A = 10-2" M^i 


t, < il,2 


-0.10 


-0.11 


0.881 


0.880 


0.23 


0.29 


M = 0c = 0.01 Mpi, M = 10=*-** Mpi, A = 10-2° M*i 


t, < il.2 


-0.17 


-0.18 


0.946 


0.955 


1.2 


1.36 


M = 0c = 0.001 Mpi,^ = 10^ Mpi, A = 10-^" M^i 


t. < il,2 


-0.14 


0.15 


0.78 


0.77 


2.9 X 10-* 


4.6 X 10-* 


M = (j}, = 0.001 Mpi, ^ = 10** Mpi, A = 10-3° M*i 


i. < ii,2 


-0.23 


-0.24 


0.77 


0.77 


0.021 


0.035 


M = 0.001 Mpi, 0c = IQ-^Mpi,^ = lO^A/puA = 10-3°M;^i 


i. < ii.2 


-0.12 


-0.13 


0.82 


0.81 


6.1 X 10-3 


8.4 X 10-3 


M = 0.001 Mpi, 0c = 10-*Mpi,/i = 10**Afpi, A = 10-3°M;^i 


t. < ii,2 


-0.20 


-0.21 


0.83 


0.84 


0.15 


0.21 



TABLE II. Comparison of the power spectrum, its tilt and the non-gaussianities in the original model for various illustrative 
points in parameter space, representing the horizon exit of the pivot scale during the qualitatively different phases. We 
compare the approximate analytical results with the numerical results obtained using the SN formalism, assuming instantaneous 
reheating. 




be neglected. Using Eq. ((T^ . one then finds 



ris - 1 



32^JW|i 
M2 



'iV.. 



M^ 
8M^ 



(108) 



This formula corresponds to the one derived in Refs. [25|, 
ISOj assuming that the waterfall trajectories are effectively 
single-field. For the amplitude of the power spectrum of 
curvature perturbations, one gets 



FIG. 1. Local /nl parameter as a function of k, for the F-term 
model, with from left to right m — 10-*Afpi, m — IQ-^Mpi, 
m — 10- Mpi, m = 10- Mpi. The bold dots are the numer- 
ical results using the SN formalism. The dotted horizontal 
line corresponds to /nl = — 5/[3(A'^oxit -l- 1)] [the approximate 
result from Eq. pOll) ] and the dashed lines to Eq. (jlOlf) . which 
are respectively valid for i« > fi,2 and £* < ti,2. We assume 
for simplicity that A'oxit ~ 60. 



1. Case 1: Hubble exit m phase 2 



In the generic case where i* > ii^2, one obtains 
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Because observable scales exit the Hubble radius near the 
critical instablitiy point, one has e* ^ 1, so that it can 
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Here also the results of Refs. [2^, [SO] are recovered. 
Therefore the regime corresponding to t* > ti^2 is ef- 
fectively single field. The F-term and D-term models 
generate a red tilted power spectrum, but the spectral 
index is too low when compared to CMB observations. 
For the original model, it is in principle possible to in- 
crease its value up to rig ~ 0.94 by increasing A that is 
a free parameter, and therefore A^oxit. However, in this 
case inflation is realised at an energy scale near the limit 
imposed by the search for i3-modc polarization in the 
CMB. 

Finally, for the original model, notice that our com- 
ment regarding the specific case where the trajectories 
reach the temporal minimum before the end of infla- 
tion still applies, and ^ond can be replaced by ^2.tm in 

Eq. mm . 
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FIG. 2. Evolution of the power spectrum of curvature pertur- 
bations (black) as well as of rescaled adiabatic (dashed) and 
entropic (dotted) perturbations, generated respectively by ini- 
tial perturbations of (j) (blue and red curves) and ip (green and 
yellow curves), and the power spectrum of curvature pertur- 
bations (black curve) for a pivot scale fc« — 0.05Mpc~^ and 
F-term potential parameters m = 10""^ and k = 5 x 10~*. 
These parameters correspond to the case i, > ii,2. 
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FIG. 3. Evolution of the power spectrum of curvature pertur- 
bations (black) as well as of rescaled adiabatic (dashed) and 
entropic (dotted) perturbations, sourced respectively by ini- 
tial perturbations of (p (blue and red curves) and tp (green and 
yellow curves), and the power spectrum of curvature pertur- 
bations (black curve) for a pivot scale k, — 0.05Mpc~^ and 
F-term potential parameters m = 10"'^ and k = 3 x 10^^. 
These parameters correspond to the case i, < ii,2 



2. Case 2: Hubble exit in 



In the case where the potential parameters are tuned 
such that t* < ii 2, we find 



diV* 

dV' 
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By using these equations as well as the relations iV_^ ^ 
^(6, N ■tb'th > N ib& and A''^^ » N &&■, we obtain the 



,i/'0 



leading term for the spectral index 
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16M|^ 
M2 
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When ^, -^ ^2 it is connected to the spectral index 
for i» > ii 2. Notice that instead of increasing when 
the Hubble exit occurs deeper in phase 1, the spectral 
index first takes lower values than expected for effec- 
tively single-field trajectories. This behaviour is shown 
in FIG. [5] that gives the spectral index as a function of 
K for the F-term model and different values of M. This 
is because |^2|, and thus |^*|, first increases with k. The 
spectral index increases up to unity only when ^, -> 0. 
From Eq. (|112p only, one can conclude that it is possible 
to find a spectral index value in agreement with CMB 
observations. It turns out however, that the power spec- 
trum amplitude is strongly modified in this case by en- 
tropy perturbations and cannot fit to CMB observations, 
as explained in the following. 

The amplitude of the power spectrum is given by 



Vc^{K 
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In the case where i* ~ ^\;ii we replace V'* — '0oexp(x2), 
such that the amplitude is 
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and it connects continuously to what is found in Eq. (jl09p 
for t^ > ti_2. When the Hubble exit of the scale k^ occurs 
deeper in phase 1, we see using Eqs. (j46|) and (|49l) that 
the amplitude grows exponentially as 



■P^(A;*) X cxp 



-(-C 



(115) 



where A^i is the number of e-folds in phase 1 (between 
the critical point and the transition to phase 2) and N^ 
is the number of e-folds between the critical point and 
the horizon exit of the scale fc* (i.e. the value of the 
parameter A^* at horizon exit, provided A^* = at the 
critical point). The spectrum then may reach a maximal 
amplitude typically larger than unity. 



Vc{U ^ 0) 
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FIGs. |4] and [5] for the F-term model illustrate this ex- 
ponential growth that prevents the amplitude to be in 
agreement with CMB observations when the spectral tilt 
is in the allowed range by CMB observations. As men- 
tioned in the previous Section, the case where the tempo- 
ral minimum is reached in phase 1 is not relevant because 
of the quantum diffusion. 
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B. From the numerical integration of multi-field 
perturbations 

We now compare with the results obtained from 
the methods explained in Section |IVl By integrating 
the multi-field perturbations, one can follow their sub- 
Hubble and super-Hubble evolution throughout the wa- 
terfall phase and identify the contributions of the fields 
to the rescaled adiabatic and entropic perturbations, that 
are respectively defined as [49| 



^TTa 



(57r„ = 



tpSij) 



tp64> (j)Stp 



(117) 
(118) 



+ -02 J(l)^+1p 



Wc have plotted on FIGs. [5] and [3] the evolution of cur- 
vature, rescaled adiabatic and entropic perturbations for 
the F-term model (but we find similar behaviors for the 
D-term and the original models) for the case t» > ti^2 
and i* < ti,2, respectively. 

In the first case, the curvature perturbations freeze 
when they become super-Hubble, as expected for an ef- 
fectively single field model. Notice also that in the sub- 
Hubble regime, the adiabatic perturbations induced by 
the initial conditions for Scj) with Sip initially set to zero 
correspond to the entropy perturbations induced by S'ip 
with S(j} initially set to zero, and conversely. This is ex- 
pected because the field perturbations evolve as indepen- 
dent plane waves with identical amplitudes in the sub- 
Hubble regime, and because 



(jTTa = cos 96(j) — sin 96ip . 
6tTc — sin d6(p + cos 9Sip . 



(119) 
(120) 



where 9 is the angle between a unit vector tangential to 
the field trajectory and the cj) = direction. 

In the second case (i» < ^1,2), we observe that super- 
Hubble curvature perturbations are receive contributions 
from entropy perturbations that are generated during the 
field evolution in phase 1 and then are frozen during the 
phase 2. This corresponds to the strong enhancement 
of the power spectrum amplitude according to Eq. (jllSp 
obtained within the SN formalism. 

We have compared the analytical approximations of 
the previous Section to the numerical results for the 
power spectrum amplitude and spectral index. For the 
F-term model, we have plotted in Figs. 0] and [S] the am- 
plitude and the spectral index of the power spectrum of 
curvature perturbations as a function of the parameter 
K, for various values of M in the F-term model. For com- 
parison, we have also plotted the effective 1-field slow-roll 
predictions and the analytical approximations derived in 
the previous Section. These plots illustrate the general 
agreement we find between the different methods, and the 
important modification of the power spectrum of curva- 
ture perturbations when compared to the predictions as- 



suming an effectively single-field dynamics. As expected, 
for the D-term model we find very similar results. 

For the original model, we have compared in TA- 
BLE [11] the power spectrum amplitude and spectral in- 
dex for some parameter sets corresponding to the differ- 
ent regimes. For the spectrum amplitude, we observe a 
strong discrepancy (up to 30%) between numerical re- 
sults and the analytical approximation in the regime 
where t,, < ti2- However, it must be noticed that a 
tiny modification of iV| — iV^,j(j can affect importantly 
the spectrum amplitude in this regime. Since this quan- 
tity cannot be evaluated analytically with a very good 
precision (even more since we assume a sharp transition 
between phase 1 and phase 2, which is not exactly the 
case), such a strong discrepancy can be explained. In 
the case t* > ti^2, numerical and analytical results agree 
well. Finally, we observe that when the temporal mini- 
mum is reached in phase 2 before the end of inflation, the 
predictions from the numerical and the approximate an- 
alytical method only show marginal deviations, because 
IC2.TMI ~ M'^/4'l ^ 1 in this particular regime. 



IX. CONCLUSION 

Using two different methods - the SN formalism and 
the numerical integration of the linear multi-field pertur- 
bations - we have calculated the level of non-gaussianities 
and the power spectrum of curvature perturbations pro- 
duced in the parametric regime of a mild waterfall phase 
in models of hybrid inflation. 

We have investigated the supersymmetric F-term 
and D-term models as well as the original non- 
supersymmetric hybrid model. In order to study these 
within a unified analysis, we have introduced a common 
parametrisation for the different variants of the poten- 
tial. For the F-term and D-term models, the mild wa- 
terfall regime occurs in the small coupling limit, when 
K < M'^/M^^. For the original model this happens when 
/iAf > Afpj. We have only considered field values lower 
than Planck mass, as it is commonly imposed for an ef- 
fective field theory description to be valid. The defin- 
ing feature of the mild waterfall regime is that the last 
60 e-folds of inflation are realized after the fields pass 
the critical point, such that possible cosmological defects 
are stretched outside the observable Universe. The ob- 
servational predictions are modified, and if one assumes 
that the waterfall trajectories are effectively single-field, 
one obtains spectral index values from rig = 1 — 4/A^oxit, 
when observable scales leave the Hubble radius in the 
so-called slow-roll phase 2 where the slope of the poten- 
tial in the valley direction is dominated by the terms 
involving the auxiliary field of the waterfall (i* > ^1,2), 
and up to unity when this happens in the first phase 
(t* < ti_2), where the slope along the inflationary valley 
dominates. Prior to the present analysis, the model was 
therefore considered as possibly in agreement with CMB 
observations [2^, |30|. However, since the scalar field tra- 
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jectories are turning during the waterfall regime toward 
the minima of the potential, one might expect a large 
contribution of entropic modes to the power spectrum of 
curvature perturbations, as well as a high level of local 
non-gaussianities. The quantification of these contribu- 
tions and the resulting signatures and the resolution of 
the question of the phenomenological viability of hybrid 
inflation in the mild waterfall regime are the main results 
of the present paper. For this purpose, we have derived 
analytical approximations for the local /nl parameter 
and the power spectrum amplitude and tilt, which agree 
well with the numerical results. 

For all the models, we find that the generic regime cor- 
responding to t* > ii,2 is effectively described by single 
field dynamics, with /nl ~ — 5/(3(A^cxit + 1)) ~ 0.03, 
which is nearly independent of the potential parameters. 
For the spectrum amplitude and the spectral index, we 
confirm the results of Refs. [2J, [2^ [SOll by using both the 
SN formalism and the linear theory of multi-field per- 
turbations. The resulting spectral index is outside and 
below the bounds imposed by CMB experiments, and 
this regime is therefore strongly disfavored. 

For the particular regime where the parameters are 
tuned so that t, < ti^2, we find that the waterfall trajec- 
tories cannot be considered as effectively single field. The 
level of non-gaussianities increases, with negative values 
of /nl down to about minus unity. Notice however that in 
no case we find that the level of non-gaussianity exceeds 
the recent bounds from the Planck experiments. Regard- 
ing the power spectrum of curvature perturbations, we 
find that the entropic modes are an important source 
of the super-Hubble curvature perturbations, enhancing 
the power spectrum by several order of magnitudes, up 
to a maximal amplitude larger than unity, which is far 
from the CMB constraints and which would lead to the 
formation of primordial black holes. When scales exit 
the Hubble radius deeper in the first slow-roll phase, the 
spectral index first takes lower values than expected for 



an effectively single field model, and then increases up to 
unity. However, values for the spectral index that are in 
agreement with CMB constraints always appear in con- 
junction with a power spectrum amplitude that is much 
higher than the measured one. Contrary to what was 
thought before [23, ^M, llOl ; it is therefore impossible to 
find parameters in agreement with CMB observations. 

We leave for a future work the particular case where 
the waterfall lasts typically 1 < A^ ^ 60, for which we 
expect a modification of the slow-roll prediction of infla- 
tion along the valley. We also would like to emphasize 
that the levels of non-gaussianities are calculated in this 
paper at the end of the inflationary era. In principle, 
they can be subsequently modified during the reheating 
era. In Ref. [56], some models are found for which the 
non-gaussianities are enhanced during the reheating. But 
in the context of a tachyonic preheating phase, applying 
the SN formalism seems to be a very challenging issue. 
Another interesting perspective would be to forecast the 
constraints on hybrid models from the CMB distortions 
generated by the perturbations from the waterfall phase 
at the end of inflation. Finally, we would like to men- 
tion that in the regime where the dynamics of both the 
fields is stochastic [57[, observable predictions have not 
been derived so far. Studying this regime will neverthe- 
less require new methods since the classical dynamics is 
not valid in this case. 
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FIG. 4. Amplitude of the power spectrum of curvature perturbations from the analytic approximations based on the SN 
formalism (Hubble exit in phase 1: dashed), from the numerical integration (points), and assuming effectively single field 
trajectories (solid) as a function of the k parameter for the F-term model, for a pivot scale fe* = 0.05Mpc~^, and from left to 
right, m = 10"VlO~VlO"V 




FIG. 5. Spectral index of the power spectrum of curvature perturbations derived from the analytic approximations based on 
the 5N formalism (Hubble exit in phase 1: dashed / phase 2: short dashed), derived from the numerical integration (points), 
and assuming effectively single field trajectories (solid), as a function of the k parameter for the F-term model, for a pivot scale 
k, = 0.05Mpc"\ and from left to right, m = 10"VlO"VlO"^ 



